Bianchi II with time varying constants. Self-similar approach. 
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We study a perfect fluid Bianchi II models with time varying constants under the self-similarity 
approach. In the first of the studied model, we consider that only vary G and A. The obtained 
solution is more general that the obtained one for the classical solution since it is valid for an 
equation of state u) G (— 1, oo) while in the classical solution ui 6 ( — 1/3, 1) . Taking into account 
the current observations, we conclude that G must be a growing time function while A is a positive 
decreasing function. In the second of the studied models we consider a variable speed of light ( VSL) . 
We obtain a similar solution as in the first model arriving to the conclusions that c must be a growing 
time function if A is a positive decreasing function. 



I. INTRODUCTION. 

Since the pioneering work of Dirac ((H), who proposed, 
motivated by the occurrence of large numbers in Uni- 
verse, a theory with a time variable gravitational coupling 
constant G, cosmological models with variable G and 
nonvanishing and variable cosmological term, A, have 
been intensively investigated in the physical literature 
(see for example 

In modern cosmological theories, the cosmological con- 
stant remains a focal point of interest (see [l5|- [l8[ for 
reviews of the problem). A wide range of observations 
now compellingly suggest that the universe possesses a 
non-zero cosmological constant. Some of the recent dis- 
cussions on the cosmological constant "problem" and 
on cosmology with a time-varying cosmological constant 
point out that in the absence of any interaction with 
matter or radiation, the cosmological constant remains 
a "constant". However, in the presence of interactions 
with matter or radiation, a solution of Einstein equa- 
tions and the assumed equation of covariant conserva- 
tion of stress-energy with a time- varying A can be found. 
This entails that energy has to be conserved by a de- 
crease in the energy density of the vacuum component 
followed by a corresponding increase in the energy den- 
sity of matter or radiation Recent observations strongly 
favour a significant and a positive value of A with mag- 
nitude A(Gh/c 3 ) « 1CU 123 . These observations suggest 
on accelerating expansion of the universe, q < 0. 

Our current understanding of the physical universe is 
anchored on the analysis of expanding, isotropic and ho- 
mogeneous models with a cosmological constant, and lin- 
ear perturbations thereof. This model successfully ac- 
counts for the late time universe, as is evidenced by the 
observation of large scale cosmic microwave background 
observations. Parameter determination from the anal- 
ysis of CMB fluctuations appears to confirm this pic- 
ture. However, further analyses seem to suggests some 
inconsistency. In particular, its appears that the uni- 
verse could have a preferred direction. Followup anal- 
yses of various sets of WMAP data sets, with different 
techniques seem to lead to the same conclusion. It is 
still unclear whether or not the directional preference 



is intrinsic to the underlying model, and what implica- 
tions this has on our understanding of cosmology. For 
this reason Bianchi models are important in the study of 
anisotropies. 

Recently , the cosmological implications of a variable 
speed of light during the early evolution of the Universe 
have been considered. Varying speed of light (VSL) mod- 
elsproposed by Moffat (|19|]) and Albrecht and Magueijo 
([2fJ|), in which light was travelling faster in the early 
periods of the existence of the Universe, might solve the 
same problems as inflation. Hence they could become 
a valuable alternative explanation of the dynamics and 
evolution of our Universe and provide an explanation for 
the problem of the variation of the physical "constants" . 
Einstein's field equations (EFE) for FRW spacetime in 
the VSL theory have been solved by Barrow ([2l| and 
[22| for anisotropic models), who also obtained the rate 
of variation of the speed of light required to solve the flat- 
ness and cosmological constant problems (see J. Magueijo 
([23l]) for a review of these theories). 

This model is formulated under the strong assumption 
that a c variable (where c stands for the speed of light) 
does not introduce any corrections into the curvature ten- 
sor, furthermore, such formulation does not verify the 
covariance and the Lorentz invariance as well as the re- 
sulting field equations do not verify the Bianchi identities 
either (see Bassett et al [24j). 

Nevertheless, some authors (T. Harko and M. K. Mak 
P.P. Avelino and C.J.A.P. Martins [H and H. Sho- 
jaie et al [13|) have proposed a new generalization of Gen- 
eral Relativity which also allows arbitrary changes in the 
speed of light, c, and the gravitational constant, G, but 
in such a way that variations in the speed of light intro- 
duces corrections to the curvature tensor in the Einstein 
equations in the cosmological frame. This new formula- 
tion is both covariant and Lorentz invariant and as we 
will show the resulting field equations (FE) verify the 
Bianchi identities. 

The purpose of this paper is to study a perfect fluid 
Bianchi II model with time varying constants through 
the self-similarity (SS) hypothesis (see foe example [341 ]- 
[42|). The study of SS models is quite important since a 
large class of orthogonal spatially homogeneous models 
are asymptotically self-similar at the initial singularity 
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and are approximated by exact perfect fluid or vacuum 
self-similar power law models. Exact self-similar power- 
law models can also approximate general Bianchi models 
at intermediate stages of their evolution. This last point 
is of particular importance in relating Bianchi models to 
the real Universe. At the same time, self-similar solu- 
tions can describe the behaviour of Bianchi models at 
late times i.e. as t — > 00. 

We would like to emphasize that in this work we are 
more interesting in mathematical respects (exact solu- 
tions) than in studying strong physical consequences. 
Nevertheless we consider some observational data in or- 
der to rule out some the obtained solutions. 

Therefore the paper is organized as follows. In the 
second section, we begin by defining the metric and cal- 
culating its Killing vectors as well as their algebra. We 
also calculate the four velocity and the quantities de- 
rived from it, i.e. the expansion, Hubble parameter etc... 
In section three, we review some basic ideas about self- 
similarity. Once we have clarified some concepts then 
we calculate the homothetic vector field as well as the 
constrains for the scale factors. In section fourth we re- 
view the classical solution i.e. a perfect fluid model with 
"constant" constants. We show that the obtained solu- 
tion belongs to the LRS Bianchi type II by calculating 
the fourth Killing vector field. We also calculate all the 
curvature invariants showing that the obtained solution 
is singular. These results will be valid for the rest of the 
studied models. We end this section showing that there 
is not solution for the vacuum model. In section five, 
we study a perfect fluid model with variable G and A. 
We compare the solution with the obtained one in the 
above section showing that we have been able to relax 
the constrains obtained for the classical model. We rule 
out some of the obtained solution taking into account the 
current observations which suggest us that A > and 
q < 0. In section six, we study a perfect fluid model with 
variable speed of light (VSL). We out line the field equa- 
tions taking into account the effects of a c-var into the 
curvature tensor. This point is essential in this approach 
as we have showed in a previous paper ([HI). We calcu- 
late the homothetic vector field for the new metric and 
we show the new constrains for the scale factors. If we 
assume a particular solution for c(t) then we show that 
it may be a decreasing time function. Only taking into 
account observational restriction we are able to rule out 
some of the obtained solutions as in the above section. 
We end summarizing the conclusions in the last section. 



II. THE METRIC. 



tively, the covariant derivative being with respect to the 
Levi-Civita connection on M derived from g. The asso- 
ciated Ricci and stress-energy tensors will be denoted in 
component form by Rij{= R c jcd) and Tjj respectively. 

A Bianchi II space-time is a spatially homogeneous 
space-time which admits an abelian group of isometries 
G3, acting on spacelike hypersurfaces, generated by the 
spacelike KVs, 



Zi=Kd x , &=d v , & = -Kyd x + d z , (1) 
and their algebra is: 
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In synchronous co-ordinates the metric is: 

ds 2 = -c 2 dt 2 + a 2 (t)dx 2 + (b 2 (t) + K 2 z 2 a 2 {t)) dy 2 + 



+ 2Ka 2 (t)zdxdy + d 2 (t)dz 2 



(2) 



where the metric functions a(t), b(t), d(t) are functions of 
the time co-ordinate only and K £ K. The introduction 
of this constant is essential since if we set K = 1, as it 
is the usual way, then there is not SS solution for the 
outlined field equations. In this paper we are interested 
only in Bianchi II space-times, hence all metric functions 
are assumed to be different and the dimension of the 
group of isometries acting on the spacelike hypersurfaces 
is three. 

Once we have defined the metric and we know which 
are its killing vectors, then we calculate the four velocity. 
It must verify, C^Ui = 0, so we may define the four 
velocity as follows: 



"' = ( ^0,0,0 i 



(3) 



in such a way that it is verified, g(u l , u l ) = —1. 
From the definition of the 4- velocity we find that: 



u 



1 



c \ a b d 
H = Hi + H2 + H3 , 

a 2 = -L {H I + H 2 + H 2 -H 1 H 2 - H X H 3 - H 2 H 3 ) 



Throughout the paper M will denote the usual 
smooth (connected, Hausdorff, 4-dimensional) spacetime 
manifold with smooth Lorentz metric g of signature 
(—,+,+,+) (see for example [33[ ). Thus M is paracom- 
pact. A comma, semi-colon and the symbol C denote 
the usual partial, covariant and Lie derivative, respec- 



III. SELF-SIMILARITY. 

In general relativity, the term self-similarity can be 
used in two ways. One is for the properties of space- 
times, the other is for the properties of matter fields. 
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These are not equivalent in general. The self-similarity in 
general relativity was defined for the first time by Cahill 
and Taub (see [23|, and for general reviews [? ]- 40]). 
Self-similarity is defined by the existence of a homoth- 
etic vector V in the spacetime, which satisfies 



(4) 



where gij is the metric tensor, Ly denotes Lie differenti- 
ation along V and a is a constant. This is a special type 
of conformal Killing vectors. This self-similarity is called 
homothety. If a =/= 0, then it can be set to be unity by a 
constant rescaling of V. If a = 0, i.e. CyQij — 0, then V 
is a Killing vector. 

Homothety is a purely geometric property of space- 
time so that the physical quantity does not necessarily 
exhibit self-similarity such as CyZ = dZ, where d is a 
constant and Z is, for example, the pressure, the energy 
density and so on. From equation (fj| it follows that 
CyR 1 jki = 0, and hence CyRij = 0,and CyGij = 0. A 
vector field V that satisfies the above equations is called 
a curvature collineation, a Ricci collineation and a mat- 
ter collincation, respectively. It is noted that such equa- 
tions do not necessarily mean that V is a homothetic vec- 
tor. We consider the Einstein equations Gij — 8irGTij, 
where is the energy-momentum tensor. If the space- 
time is homothetic, the energy-momentum tensor of the 
matter fields must satisfy CyTij — 0. For a perfect 
fluid case, the energy-momentum tensor takes the form 
Tij = (p + p)uiUj + pgij, where p and p are the pressure 
and the energy density, respectively. Then, equations (|4|) 
result in 



Cvu 1 = 



Cyp — —2ap, Cvp — ~2ap. (5) 



As shown above, for a perfect fluid, the self-similarity 
of the spacetime and that of the physical quantity co- 
incide. However, this fact does not necessarily hold for 
more general matter fields. Thus the self-similar vari- 
ables can be determined from dimensional considerations 
in the case of homothety. Therefore, we can conclude ho- 
mothety as the general relativistic analogue of complete 
similarity. 

From the constraints (O, we can show that if we con- 
sider the barotropic equation of state, i.e., p = f(p), then 
the equation of state must have the form p = top, where u> 
is a constant. In this paper, we would like to try to show 
how taking into account this class of hypothesis one is 
able to find exact solutions to the field equations within 
the framework of the time varying constants. 

From equation (f5]),we find the following homothetic 
vector field for the BII metric @: 



v = td t + \1 - t-j xd x + [l - t-j yd y + [l - t-j zd z , 

(6) 

with the following constrains for the scale factors: 

a(t) = a t ai , b(t) = b t a2 , d(t) = d Q t a3 , (7) 



with a%, <22, ct3 £ M, in such a way, that the constants 
must verify the following restriction: 

a-2 + a 3 - ai = 1. (8) 

therefore the resulting homothetic vector field is: 

V = td t + (1 - at) xd x + (1 - 03) yd y + (1 - a 3 ) zd z , (9) 

which is the same as the found one in the Bianchi I model 
(see for example (43|). 

IV. THE PERFECT FLUID MODEL. THE 
CLASSICAL SOLUTION. 

Taking into account the field equations (FE) 



(10) 



2 — e* 

where, in this section, we shall consider that A vanish, 
and the energy- momentum tensor, T^, is defined as fol- 
lows: 



Tij = (p+p)utUj + pg^, 



(11) 



with the usual equation of state p — u>p, to £ R, as we 
have discussed above. 

The resulting FE for the metric © with a perfect fluid 
matter model (fTTj) are: 
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-up, 



(12) 
(13) 
(14) 



— (15) 



, a' b' d! \ 
(f+p(l + w) + - + =0. 



(16) 



Now, if we take into account the obtained SS restric- 
tions for the scale factors i.e. 

a(t) = a t ai , b(t) = b t a2 , d(t) = d t as , 

where 01,02,03 £ K, must satisfy 02 + 03 — 01 = 1, we 
get from eq. (Q2 



P = Pot-\ (17) 

where 7 = a (1 + u>) = (at + 02 + a 3 ) (1 + u>) . Therefore 
the system to solve is the following one: 

3A- 2 



a 2 (a 2 - 1) + a 3 (a 3 - 1) + a 2 a 3 



4 

K 2 



= -Auj, (18) 



a 3 (a 3 - 1) + aia 3 + at (at - 1) + — = -Auj, (19) 



02(02 - 1) + ai«2 + ai(ai - 1) + 



4 



at + a 2 + a 3 



-Auj, (20) 
2 (21) 



1 + oj' 
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where we have set the constant A as follows 

A = a x a 2 + aia 3 + a 3 a 2 — , 

and 02 + as — ax = 1- 

The obtained solution is the following one: 



1 



ax 



2(w + l)' 



a-2 = «3 



therefore the constrain a 2 
2a 2 — cl\ = 1, and 



«3 



_ w + 3 
" 4 + 1) ' 

ai = 1 collapses to 



A" 



j _ 1 + 2w - 3w 2 



4(w + l) z 



1 



-, 1 



soa 2 e(±,l), ai £ (0,1), K 2 G (0, 1/4) , and A = 
2aia2 +a^-4" = ■ Taking into account all these 

results the homothetic vector field collapses to 

V = td t + (1 - oi) x9 x + (1 - a 2 ) (yd y + zd z ) , 

this is the already solution obtained by Hsu et al (Ell). 
This solution has been already obtained by Collins (|44j) 

Remark 1 As is it observed, we may find the relation 
a 2 + a 3 — a i = 1 from dimensional considerations, only 
looking at the FE. Note that the quantity 



a 2 c 2 
b 2 d 2 



2+2ax 



aft- 



b\t 2ai >dlt 2a * 



t 



2a\ — 2a 2 — 2a3 



after substitution, it must have dimensions ofT 2 as the 
rest of the factors. Therefore we find the relation a 2 + 
03 — ai = 1. 

Therefore we have obtained a LRS BII model, note 
that d = b, and hence the metric collapses to this one: 

ds 2 = -c 2 dt 2 + a 2 (t)dx 2 + (b 2 (t) + K 2 z 2 a 2 (t)) dy 2 + 
+ 2Ka 2 (t)zdxdy + b 2 {t)dz 2 , (22) 

which admits the three Killings vectors ([I]) and this new 
one 



-K 



Hence their algebra is: 



■L 
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d x + zd y - yd z 



(23) 
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i.e. Cgg = 1, 



C24 



T, and C| 4 = 1 



With the obtained results we can see that 
t V w + 1 / * 

while the shear behaves as: 

2 (a 2 - l) 2 / 3cj + 1 \ 2 1 



3c 2 t 2 V 4 (^ + 1 )/ 3c 2 £ 2 ' 

As it is observed, these quantities fit perfectly with the 
current observations by High -Z Supernova Team and 
Supernova Cosmological Project see for example: Gar- 
navich et ai, 1998 Perlmutter et ai, 1997, 1998, 

1999 [111; Riess et ai, 1998 [U; Schmidt et ai, 1998 
EH- 



A. Curvature behaviour. 

With all these results, we find the following behaviour 
for the curvature invariants (see for example |45l| - [jj| ) ■ 
Rice Scalar, I , yields 



2 / A 2 c 2 
/ o = ^72 (lla!-10a 2 + 2- — 



(24) 



while Krestchmann scalar, I\ := RijkiR 1 ^ 1 , yields: 

h = — j-j (432a 2 - 960a^ + 896a 2 , - 384a 2 + 64 
+A 2 c 2 (HA 2 c 2 - 40a 2 , + 80a 2 - 48)) . 
The full contraction of the Ricci tensor, I 2 := RijE?° , 



is: 



h = 



(528a 2 - 1024^ + 768a 2 2 - 256a 2 + 32 



4c 4 t 4 

K 2 c 2 (3A 2 c 2 - 16a 2 + 8)) , 

this means that the model is singular. 

The non-zero components of the Weyl tensor. The 
following components of the Weyl tensor run to ±00 when 
t 0, 

n ~ C ~ +4(a2-l) C ^ x2(ai— 1) 
^txtx ~ ^txty ~ I , ^tztz ~ I j 



and 



these others depend on the value of a 2 , 

/~> ~r i ~ ~ +4a 2 — 3 

^txyz ~ ^tyxz ~ ^tyyz ~ ^tzxy ~ ^ ) 



/"< ~ ~ /° ~ j-2(3a 2 — 2) 

^121/2 ~ ^3:2:3:2: ~ ^xyxy ~ j 

,o ^ ,2(2a 2 -l) , 2,2(3a 2 -2) 

^yzyz ~ h ~r 6 b 
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i.e. they may run to zero as well as to ±00 when t — > 0, 
depending on the value of a 2 . 

The Weyl scalar, I 3 = C abcd C abcd = h - 27 2 + \l 2 , 
(this definition is only valid when n = 4) 

4 



([(03 - f ) (-2 (o 2 - 1) - 3Xc) + if 2 c 21 



3c 4 t 4 



[(a 2 - 1) (-2 (oa - 1) + 37sTc) + if 2 c 2 ]) 



The electric part scalar 74 = EijE^ , (see W.C. Lim et 
al E|) 



i 



6c 4 t 4 



-2(a 2 -l) 2 + /f 2 c 2 ) , 



while the magnetic part scalar I5 = H^H^, yields 

Therefore the Weyl parameter W 2 (see W.C. Lim et 
al [H) 

therefore 

W2 = 36ct ( ( ° 2 1)2 + K2c2 T ( (2 ° 2 " 2)2 + ^ 2 

note that the value of W 2 is reall y sm all. 
The gravitational entropy (see [4(|-[47j]) 



p2= h = h - 2/ 2 + \l 2 _ h , 1 IJ 



2 = const. 



h h 3/ 2 

since the spacetime is SS (see [5(| and [43[ for a discus- 
sion) . 



Therefore, we conclude that there is no physical so- 
lution for the vacuum model. Note that when we set 
K = 0, the metric is reduced to Bianchi I one. In the 
same way, it is also verified the relation 02+03 — 01 = I. 



V. PERFECT FLUID MODEL WITH G AND 
LAMBDA VARIABLE. 

The resulting FE for the metric © with a perfect fluid 
matter model (jlip and with G and A time-varying are: 



a' b' 


a' d! 


d! b' 


K 2 


2 2 
a c 


~T + 

a 


a d 


1 d b 


4 


b 2 d 2 


b" 


d" 


d' b' 


ZK 2 


a 2 c 2 




' ~~d 4 


~d~b ~ 


4 


b 2 d 2 


d" 


a" 
+ 


a' d' 


K 2 


a 2 c 2 


d 


a 


a a 


+ 4 


b 2 d 2 


b" 


a" 


a' b' 


K 2 


a 2 c 2 


b 


H 

a 


a b 


+ 4 


b 2 d 2 



8ttG 



p + Ac 2 , (29) 



8ttG 



-up + Ac 2 , 
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(32) 
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(34) 



Now, we shall take into account the obtained SS re- 
strictions for the scale factors i.e. 



B. Vacuum model. 

In this section we shall show that there is no self-similar 
solution for the vacuum model. In this case, after sub- 
stitution, the resulting FE are: 
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(28) 



so, by solving this system, we only obtain the following 
unphysical solutions: 

«i = o 2 — a 3 = K = 0, 

01 = 1, a 2 = a 3 = K = 0, 



Ol 



1 1 ^2 4 

02 = a 3 = -, K = --, 



3' " 3 

01 = — -, o 2 = a 3 = -, K = 



a(t) = a a t ai , b(t) = b a t a2 , d{t) = d Q t a \ 

where 01, a 2 , 03 G K, in such a way that they must verify 
02 + o 3 — Ol = 1. 

From eq. (f3"3"]l we get 



(35) 



P = Pot 7 , 

where 7 = (u + 1) a and a = (a\ + 02 + 03). 
From eq. ([29)) we obtain: 



BttG 



-Po 



(36) 



where A = 0102 + 0103 + 0203 , . 

Now, taking into account eq. ([34]) and eq . ([36]) . algebra 
brings us to obtain 



G = Got 7 " 2 , 



where Gn = -. — i A , t\ ■ 

While the cosmological "constant" behaves as: 



(37) 



A-4ll- 
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= A r 



(38) 
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With all these result we find that the system to solve 
is the following one: 



012(02 - 1) + 03 (a 3 - 1) + a 2 a 3 - 



3K 2 



= A 



a-2 



(39) 



s K 2 a (ol-2 
«3 («3 - 1) + ai«3 + ai (ai - 1) + — = A 



(40) 



K 2 



a 2 (a 2 - 1) + aia 2 + ai(a ± - 1) + — — = ^4 

4 \ a 



a-2 



(41) 



whose solution is: 



ai = 2a 2 — 1, a 2 



03, 



-4«2 + 6a 2 - 2, (42) 



therefore this solution has only sense if a 2 G (|, l) . Note 
that for these values, a\ > 0, a\ = 2a 2 — 1 G (0, 1) . This 
result is valid for all equation of state. We would like 
to point out that this solution is more general that the 
obtained one in the perfect fluid case with "constant" 
constants, since it is valid for lo G ( — 1, 00) , although we 
are only interested in lo G (— 1, a] , a > 1. It is important 
to emphasize that w ^ -1. 

With this solution we may see that all the quantities 
depend on two variables, a 2 and lo. Newton gravitational 
constant behaves as follows: 



G 



c 2 A 



4-7T pa (lo + I) a 



f- 2 = Got 9 , Go > 0, (43) 



with g = 

2a~ K 



2 G (-2,4) , and A = ha\ - 
G (|,3) . So depending of 



(W + I)(4fl3-1) 

i 2 — =j- = 6a 2 — \a 2 + g 
the different combinations between a 2 and lo, G may be 
a growing or a decreasing time function. 

With regard to the behaviour of the cosmological con- 
stant, as we can see, it is a decreasing time function, as 
it is expected. Only rest to know which is the sign of Ao. 
We may express its behaviour in the following tables 
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/A 
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/I 



this means, for example, that when we fix a 2 \ 1/2 i.e. 
that a 2 tends to 1/2, for all value of lo between -1 and 
1, A < 0, A = 0, if uo = 1 and only A > when 
lo > 1. In the say way we may say that if we fix a 2 /* 1, 
then G (-1, -1/2) , we find that A < 0. A = 0, if 
lo = —1/2 and only it is found Ao > when lo > —1/2. 
So may only say that the sign of Ao, depends on the 
parameters a 2 and lo. If we take into account the current 
observations which suggest that Ao > 0, then we may 
rule out the values which make Ao < 0. 

As it is observed the behaviour of G and the sign of A , 
are related. If G is growing then Ao > 0, if G — const. 



then Ao = 0, i.e. the model collapses to the standard 
one studied above, and if G is decreasing then A < 0. 
Therefore if we consider the recent observations which 
suggest us that Ao > 0, we must rule out the other cases, 
concluding that G is a growing time function and A is a 
positive decreasing function. 
The energy-density behaves as 



P = Pot 7 , 



with 



7 6(0,6) 



(44) 



so it is always a time decreasing function if a; > —1, and 
to end, we find the following behaviour for the Hubble, 
deceleration and shear parameters: 



H = 



^a 2 - 1 



4a 2 



4a 2 



<0, 



(°2 - 1) 

3c 2 t 2 



0. 



(45) 



As in the above model, the metric of this one collapses 
to a LRSBII type 



VI. PERFECT FLUID MODEL WITH VSL. 

We start by defining the new metric as follows: 
ds 2 = ~c(t) 2 dt 2 + a 2 (t)dx 2 + (b 2 (t) + K 2 z 2 a 2 {t)) dy 2 + 



+ 2Ka 2 (t)zdxdy + d 2 (t)dz 2 



(46) 



note that we have simply replaced c by c(t), so we shall 
consider the four velocity: 

u= ^-L,0,0,0j, / (47) 

The time derivatives of G, c and A are related by the 
Bianchi identities 



1 



-Rgi. 



•j 



1 -^r T ii ~ A.'/,., 



(48) 



•2"-'" J V (■ 
in our case we obtain 

, + , a+u)H+ *g + ,(!_^)_c (49) 

where H = — + 7- + 4-, but if we take into account 

a b a ' 

the condition, Tf- — 0, it is obtained the following set of 
equations: 



p'+p(l + Lo)H = 0, 

8nGp + G c ~ ' 
Therefore, the resulting FE are: 



(50) 
(51) 



7 



~b ' 

<i 
T 



a' b' a' d' 
a b ad 



dT_ d/_tf 
d + d b 
a/_d/_ a" 
ad a 

a' 
a 



a b 



c \b 



a 



<2V_ 
' ~d~b 
df_ 
d 
d[ 

~i 
v_ 

1) 



K 2 a 2 c 2 

~ ~¥d 2 

3K 2 a 2 c 2 
4 b 2 d 2 
K 2 a 2 c 2 



a (I i 4 b 2 d 2 
a is \ K 2 a 2 c 2 
a h ) 4 b 2 d 2 



. a ; b' d' 



A'c 4 G' d 

1 4— 

8ttG P G c 



8ttG 
c z 

8vrG 

—P- 

c A 

8nG 
8-kG 



0. 
0. 



Ac 2 



Ac 2 , 



Ac 2 , 



Ac 2 , 



(52) 
(53) 
(54) 
(55) 
(56) 
(57) 



Note that we have taken into account the effects of a 
c— var into the curvature tensor (see [28j for a discussion 
about the issue). 

In this model, the homothetic vector field is: 



X 



dt 



Jcdt 
f cdt 

1 - ' ^ 



with the following restrictions 



[cdt 
c(t) 



f cdt 



(58) 



a = ao 



cdt 



b 



cdt 



d 



cdt 
(59) 

with (cn) i=l G K, in such a way that they must satisfy 
the relation a.2 + a 3 — ct\ = 1. Since we expect to get 
a growing scale factors, we shall consider that (/ cdt) 
must be a positive growing time function. Note that c 
only needs to be integrable but it may be decreasing as 
we shall show below. 

Remark 2 As it is observed, ifc = const., we regain the 
usual homothetic vector field, i.e. 

X = td t + (l- tHx) xd x + (1 - tH 2 ) yd y + (1 - tH 3 ) zd z , 

while the scale factors behave as 

a = a (t) ai , b = b (t) a \ d = d {t) a \ (61) 
with tt2 + «3 — ctx = 1, as in the above studied cases. 
By defining 



H, 



'Jcdt 



H = a 



Jcdt 1 



(62) 



where a — Y^i=i a i-> we nn d, from eq. (|56[) . the behavior 

cdt] , (63) 



of the energy density i.e 
P = Po 



where 7 = (1 + u>)a. 

In the same way it is easily calculated the shear 



1 



3c 2 (t) 



(H 2 + Hi + H 2 - H X H 2 - HiH 3 - H 2 H 3 ) 



\»=1 i& j 



cdt 



(64) 



As it is observed all the quantities depend on J c(f)dt. 
Now only rest to calculate G and A. 
From eqs. (1521 and |6"3"|) we get: 



8ttG 



-po 



Ac 2 



(65) 



where we have written, for simplicity, J c instead of J cdt, 
and we have set A — a±a 2 + cticx.3 + a 2 otz — j-, therefore 



A' 



2 Ac 8tt Pq G 



(fey cHj c y 

Now, taking into account eq. (15"7| , we get that 



G'd c 
_4-- 7t - 
G c J c 



(66) 



8nGp 



2Ac 



8irpoG 



G' d c 
— -4- - 7t - 
G c J c 



G' d 

G- 4 c^ ' 



and hence we obtain 



G = 



.4 



4-7rpo7 



7-2 



(67) 



(68) 



and in this way we find that the cosmological constant 
behaves as 



A = A\ 1- - 

7 



(69) 
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As we can see, from eqs. (1681 and 
verified the following relationships 



we have that are 



Gp 



A 



= const. 



(70) 



Now, we will try to find the value of the constants 
(afj)j_j . Taking into account the field eqs. (|52M55[) we 



find that, obviously eq. (|52[) vanish, but from eqs. 
155]) we get 



a 2 (a 2 - 1) + a 3 (a 3 - 1) + a 2 a 3 



3K 2 



K 2 



«3 (a.3 - 1) + aia 3 + ax (ai - 1) + — 



K 2 

a 2 {a 2 - 1) + aia 2 + ai(ai - 1) + — 



.4 



= A 



a - 2 



(71) 



a - 2 



(72) 



a - 2 



a 



(73) 



finding therefore the same solution as in the above model, 
i.e. 



a i 



2a 2 — 1, a 2 = a 2 , AT 



-4a^ + 6a 2 -2, (74) 



this solution has only sense if a 2 6 (5,1), note that for 
these values, a\ > 0. 

We would like to emphasize that to out line the system 
(|7H73j) it is essential to take into account the effects of 
a c-var into the curvature tensor. If we do not consider 
such effects then the system depends on J c and therefore 
it is not algebraic (see [28| for a detailed discussion). 

Therefore the behaviour of the main quantities is the 
following one. The energy density behaves as 



P = Po 



• (t) dt 



(75) 



where 7 = {uj + 1) a, and a = 2a 2 + a±, while Hubble 
and the deceleration parameters behave as 



H = a- 



Jcdt' 
and the shear yields 



Q I C C 



1, 



(76) 



a 2 = l -{ a2 -lf (J cdt 



On the other hand, G behaves as, G = 
G0C 4 (/ cdty , where Go > 0. Hence, its behaviour will 
depend on the different values of 7. If 7 E (0, 2) we find 
that G is a decreasing time function. If 7 = 2, G = const. 
and all the model collapses to the standard one studied 
in the above sections (G = const., c = const, and A = 0). 
While if 7 E (2, 6) then G is a growing time function. 



The cosmological constant behaves as 
2 



A = A1 



7 



-2 , . N -2 

cdt ) = A ( / cdt 



where the sign of Ao depends on the different values of 

7, 

' <0 if 7 6(0,2) 
Ao = { = if 7 = 2 

>0 if 7 6(2,6) 

Since our model is formally self-similar, then ([37l|-(36l|) 
have shown, that all the quantities must follow a power 
law, so, we may assume that for example, c takes the 
following form: c{t) = cot e , with e E K. Hence, the scale 
factors behaves as: 



a = a t ai( - e+1 \ b = M Q2(c+1) = d, 



(77) 



in such a way that they behave as growing time function 
if e E (— l,oo) , ai > 0. This means that c may be a 
decreasing time function. The nomothetic vector field 
collapses to this one 



V 



1 



d t + (1 - a 2 ) (2xd x + yd y + zd z ) (78) 



The Hubble parameter behaves as: H = (4a 2 — 
l)(e + l)t _1 = at^ 1 , so we find again the restriction 
for e, e E (— 1, 00) . In the same way, we find from eq. 
A = Ao (/ c) , that, since A must be a decreasing time 
function, this is only possible if e E (— 1, 00) , the special 
case, e = —1, is forbidden, note that f cdt — -^jt e+1 > 
0, Ve E (-l,oo). 

So, we find that 



P ■ 



£-70+1) Q ~ £7(«+l)+ 



2e-2 



-2(e+l) 



(79) 



with e E (—1, 00) , and 7 = (ui + 1) a = (u> + 1) (4a 2 — 1). 
Therefore the energy density is always a decreasing time 
function if co E (—1, 1] and e E ( — 1, 00) . As above, the 
behaviour of G depends on (a 2 , to, e) , hence it may be a 
growing function as well as a decreasing one. With regard 
to the cosmological constant, the most important thing 
is to know its sign. Recalculating all the computations 

we find that A = Ac^ 2 ( 1 — 1^, where 7 = 7 (e + 1) = 

(uj + 1) (4a 2 — 1) (e + 1), in such a way that this sign is 
affected by the perturbing parameter (e + 1) . We have 
set A — (5a 2 — 2a 2 ) (e + l) 2 — remember that, a% = 
2a 2 - 1. 

To clarify this results we may fix u> — 1, a 2 — 2/3 and 
e = ±1/2, finding 



e 


P 


G 


A 


-1/2 


f -5/3 


£-4/3 


< 


1/2 


t~ 5 


i 4 


> 
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Hence, if we take into account the current observation, 
then we must rule out those values which make A < 0. 
Therefore we conclude that c and G are growing time 
functions while A is a positive time decreasing time func- 
tion. 

With regard to the deceleration parameter we find that 
it has the following behavior: 

<7=-l + 4, (80) 
a 

where a = (4a2 — 1) (e + 1) , indicating as that it is quite 
unlikely that e — > — 1. We find that q < 0, if a > 1. 

VII. CONCLUSIONS. 

We have studied several perfect fluid Bianchi II mod- 
els under the self-similarity hypothesis. We have started 
our study reviewing the classical solution and showing 
that it belongs to the LRS BII kind by calculating the 
fourth Killing vector. The self-similar solution brings us 
to get a power law solution for the scale factors i.e. they 
behave as a = aot ai , b — b t a2 = d, in such a way the 
the constants aj must satisfy the relation 2a2 — a\ = 1. 
The only drawback is that this solution is only valid for 
an equation of state u 6 (—1/3, 1). We have also shown 
that this solution fits perfectly with the current observa- 
tions for H and q. With regard to its curvature behaviour 
we may conclude that the model is singular. The Weyl 
parameter is quite small and the gravitational entropy is 
constant as it is expected for a self-similar space-time. 
This happens because the definition of gravitational en- 
tropy does not work in this kind of space-times. We have 
finished our study of the classical model showing that 
there is not solution for the vacuum model. 



With regard to the perfect fluid model with G and A 
variable we have arrived to the conclusion that its solu- 
tion is quite similar to the obtained one for the classical 
model. But we have been able to enlarge the set of values 
for the equation of state, in this case is w £ (— 1, oo) . As 
we have mentioned above the solution is similar to the 
classical model so it also belongs to the LRS BII kind. 
With regard to the behaviour of G and A we have showed 
that G is a growing time function while A is a positive 
decreasing time function. We have ruled out the rest of 
the solution taking into account the recent observation 
which suggest us that Ao > and q < 0. 

In the VSL model we have arrived to similar conclu- 
sions as in the above model i.e. the solution is valid 
for li £ (—1, oo). We have started outlining the FE but 
taking into account the effects of a c— var into the cur- 
vature tensor. In this approach is essential to take into 
account such effects, in other way, it is impossible to get, 
after substitution, and algebraic system of equations. We 
have also calculated the homothetic vector field for the 
new metric as well as the constrains for the scale factors. 
In this case all the quantities depend on J c(t)dt. By solv- 
ing the associated algebraic system we have obtained the 
same solution as in the above studied case. Hence this 
solution belongs to the LRS BII type. In a generic way 
we have discussed the behaviour of the main quantities 
ruling out solutions like A < 0. If we assume a power 
law for c(t), c = t e say, then we arrive to the conclu- 
sion that it may be a decreasing time function. The only 
restriction is J c(t)dt must be a growing time function, 
hence c must be integrable but it is allowed that it may 
be decreasing. But if we consider the possibility of a c de- 
creasing then Ao < 0. Therefore we have concluded that 
G and c are growing time functions while A is a positive 
decreasing time function. 
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